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Abstract. We propose helium-4 spallation processes induced by long-lived stau in 
' supersymmetric standard models, and investigate an impact of the processes on light elements 

abundances. We show that, as long as the phase space of helium-4 spallation processes is open, 
they are more important than stau-catalyzed fusion and hence constrain the stau property. This 
talk is based on works pQ. 
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Q\ 1 1. Introduction 

Long-lived charged massive particles (CHAMPs) will play interesting roles in the Big-Bang 
Nucleosynthesis (BBN). The light nuclei will interact not only with the CHAMPs during the 
BBN processes El El ED El El E], but also with the decay products of the CHAMPs in the 
post-BBN era [HIE]. The standard BBN will thus be altered, and so is the abundance of the 
^ " light elements at the present time. One can thus constrain the models beyond the Standard 

Model by evaluating their prediction on the light elements abundance and comparing it with the 
current observations. We can then give stringent predictions for the forthcoming experiments 
and observations according to these constraints. 

The Standard Model extended with supersymmetry (SUSY) is one of the models that can 
accommodate such long-lived CHAMPs. With the i?-parity conservation, the lightest SUSY 
particle (LSP) is stable and become a cold dark matter. Interestingly, it can offer a long-lived 
CHAMPs if the LSP is the bino-like neutralino Xi- Coannihilation mechanism is required to 
account for the dark matter abundance in this case |10| . where the LSP and the next-lightest 
SUSY particle (NLSP) are almost degenerate in mass. Staus, denoted by f and a possible 
candidate of the NLSP, can acquire a long lifetime when the mass difference with the LSP is 
less than the mass of tau leptons. This is due to the phase space suppression of the final state 
that necessarily consists of three particles or more. Noting that such long-lived staus will be 
copious during the BBN [Til 02] , we have shown in [31 El [7j that their presence indeed alters the 
prediction of the standard BBN and possibly solve the discrepancy of the lithium abundance in 
the Universe through the internal conversion reactions. 



In this talk, we propose new reactions 



(f 4 He) 
(f 4 He) 
(f 4 He) 



x\ + v T + t + n, 
X? + v t + d + n + n, 
Xi + Vt + P + n + n + n, 



(la) 
(lb) 
(lc) 



in which (f 4 He) represents a bound state of a stau and He nucleus. Reaction (pTJ) is essentially 
a spallation of the 4 He nucleus, producing a triton t, a deuteron d, and neutrons n. Presence of 
such spallation processes has been ignored so far due to the naive expectation that the rate of 
the stau-catalyzed fusion [2] 

(2) 



(f 4 He) + d^ f + b Li 



is much larger than that of 



is larger than the reaction ([I]). Indeed, the cross section of Eq. 
4 He + d — > 6 Li + 7 by (6 - 7) orders of magnitude [13]. 

We point out that this expectation is indeed naive; the reaction Eq. (prj is more effective 
than Eq. ([2]) as long as the spallation processes are kinematically allowed. The former reaction 
rapidly occurs due to the large overlap of their wave functions in a bound state. On the other 
hand, the latter proceeds slowly since it requires an external deuteron which is sparse at the 
BBN era. The overproduction of t and d is more problematic than that of 6 Li. This puts new 
constraints on the parameters of the minimal supersymmetric standard model (MSSM). Note 
that there is no reaction corresponding to Eq. ([1]) in the gravitino LSP scenario |14t 1 15 j . 

2. Spallation of helium 4 

The 4 He spallation processes of Eq. ([1]) is described by the 
Lagrangian 

C = f *x?(3lPl + 9r P r )t + y/2G F v T -fPiT J M + h.c, 

(3) 

where Gp = 1.166 x 10 -5 GeV -2 is the Fermi coupling 
constant, -Pl(R) represents the chiral projection operator, 
and J M is the weak current. The effective coupling 
constants #l and g R are given by 
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(4) 



where g is the SU (2)l gauge coupling constant and 6*w is 
the Weinberg angle. The mass eigenstate of staus is given 
by the linear combination of tl and tr, the superpartners 
of left-handed and right-handed tau leptons, as 



f = cos y T TL + sin 9 T e 



-1-fT 



TR- 



(5) 



Here T is the left-right mixing angle of staus and 7 T is 
the CP violating phase. 




Figure 1. 4 He spallation processes. 



2.1. (f 4 He) -^Xi + Vr + t+n 

First we consider the process of Eq. (jlap . The rate of this process is expressed as 

r tn \ip\ 2 ■ ^ 

where \ip\ 2 stands for the overlap of the wave functions of the stau and the 4 He nucleus. We 
estimate the overlap by 

7T 

where Z and mn e represent the atomic number and the mass of 4 He, respectively, and a is the 
fine structure constant. We assumed that the stau is pointlike particle and is much heavier than 
4 He nucleus so that the reduced mass of the bound state is equal to the mass of 4 He nucleus 
itself. The cross section of the elementary process for this reaction is denoted by avt n and 
calculated as 



av tn = o-u((f 4 He) -)• xVrtn) 



2E f J (2tt) 3 2E u (2vr) 3 2^ (2vr) 3 (2vr) 3 (8) 
x |.M((f 4 He) ^xMf 

X (2vr)V 4) {p f + p Ue -Pu-Qt- <7n)- 

Here pi and E{ are the momentum and the energy of the particle species i, respectively. 

We briefly show the calculation of the amplitude of this process. The amplitude is 
deconstructed as 

.M((f 4 He) -> x?^tn) = (tnx?^|Ant| 4 Hef) 

= (tn|JH 4 He) (x>r\j„\f). 

The weak current J M consists of a vector current and an axial vector current as 
Ju, = Vu, + 9A^fi, where g A is the axial coupling constant. The relevant components of the 
currents in this reaction are V° and A 1 (i = 1,2,3). We take these operators as a sum of a 
single-nucleon operators as 

4 4 

V ° = E T « eiq ' r " > A * = E T a<j qr " , (10) 

a=l a=l 

where q is the momentum carried by the current, r a is the spatial coordinate of the a-th nucleon 
(a G {1, 2, 3, 4}), and r~ and a % a denote the isospin ladder operator and the spin operator of the 
a-th nucleon, respectively. Each component leads to a part of hadronic matrix element: 

(tn|y°| 4 He) = V2~M tn , 
(tn| 0A ^ + | 4 He) = V2g A M tn , 
(tn| 0A ^"| 4 He) = -V2g A M t n, 



(tn| OA ^ 3 | 4 He) = -V2g A M 



til ; 



Table 1. Input values of the matter radius R ma t for d, t, and 4 He, the magnetic radius R 
for p and n, nucleus mass mx, excess energy Ax for the nucleus X, and each reference. 



nucleus 


-Rmat(mag) [fm]/[GeV l ] 


m x [GeV] 


Ax [GeV] 






P 


0.876 / 4.439 [E] 


0.9383 [21] 


6.778 x 10" 


-3 


[22] 


n 


0.873 / 4.424 [T7] 


0.9396 [21] 


8.071 x 10" 


-3 


[22 


d 


1.966 / 9.962 p2] 


1.876 [22] 


1.314 x 10" 


-■> 


22 


t 


1.928 / 9.770 [H] 


2.809 [22] 


1.495 x 10" 


-■> 


22 


4 He 


1.49 / 7.55 [20] 


3.728 [22] 


2.425 x 10" 


-3 


[22] 



where A ± = (A 1 ± iA 2 )/y/2. Given the relevant wave functions of a 4 He nucleus, a triton, and 
a neutron (see Appendix in Ref. [I]), we obtain the hadronic matrix element as 



M 



tn 



1287T 



OHeOt 



3 (a Hc + «t) 



3/4 



jexp 


[ ^ ' 


— exp 


3an e _ 







(gt + 9n) 2 

3an e 6(a Hc + at) 



(12) 



Here q t and q n are three-momenta of the triton and the neutron, respectively, and an e and at 
are related to the mean square matter radius R m by 
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(13) 



We list in Table Q] input values of the matter radius for the numerical calculation in this article. 
The remaining part is straightforwardly calculated to be 
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Kx?^|j,|f)| 2 = 4G2| OR | 2 
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(14) 



where E v and pf, are the energy and the z-component of the momentum of the tau neutrino, 
respectively. We assumed that the stau and the neutralino are non-relativistic. This equation 
includes not only all the couplings such as Gp, gi,, and or, but also the effect of the virtual 
tau propagation in the Fig. 1. Note here that ol coupling does not contribute. This is because 
the virtual tau ought to be left-handed at the weak current, and it flips its chirality during the 
propagation since the transferred momentum is much less than its mass. 

Combining hadronic part with the other part, we obtain the squared amplitude as 



|7W((f 4 He) ^xiMn) 



2 8m^oGl\g R \ 2 



mi 



{l + 3,g\)E v \Mt 



(15) 



Integrating on the phase space of the final states, we obtain the cross section as 
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(16) 
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Figure 2. Factor I tn in Eq. f)16|) as a function 
of mass difference between the stau and the 
neutralino. Here we took mf = 350GeV, 
sin# T = 0.8, and j T = 0. 
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Figure 3. Timescale of spallation processes as 
a function of 5m and the stau-catalyzed fusion 
at the universe temperature T = 30keV |13j . 
The lifetime of free f (solid line) is also 
depicted. Here we took mf = 350GeV, sin 9 T = 
0.8, and 7t = 0. 



Here Itn is the numerically calculated factor including the information of phase space for this 
process. Analytic form of it n is shown in Ref. [T], and numerical result is depicted in Fig. 2. 
A tn , k±, and k n are defined as 

A tn = 5m + A Hc - A t - A n - Eb, 

/ / ( 17 ) 

k t = V2m t A t n , K = V2m n A tn , 

where Ax is the excess energy of the nucleus X, and E^ is the binding energy of (f 4 He) system. 

3. Comparing the rate of spallation reaction with that of stau-catalyzed fusion 

We compare the rate of the spallation and that of the stau-catalyzed fusion. We first note 
that the rate of stau-catalyzed fusion strongly depends on the temperature [T3], and we fix 
the reference temperature to be 30keV. Staus begin to form a bound state with 4 He at this 
temperature, which corresponds to cosmic time of 10 3 s. Thus the bound state is formed when 
the lifetime of staus is longer than 10 3 s. 

Figure O shows the timescale of the spallation processes as a function of 5m. The lifetime of 
free stau is plotted by a solid line. We took the reference values of mf = 350GeV, sin# T = 0.8, 
and j T = 0. The inverted rate of the stau-catalyzed fusion at the temperature of 30keV is also 
shown by the horizontal dashed line. Once a bound state is formed, as long as the phase space 
of spallation processes are open sufficiently that is 5m > 0.026GeV, those processes dominate 
over other processes. There f property is constrained to evade the over-production of d and/or 
t. For 5m < 0.026GeV, the dominant process of (f He) is stau-catalyzed fusion, since the free f 
lifetime is longer than the timescale of stau-catalyzed fusion. Thus light gray region is forbidden 
due to the over-production of 6 Li. 

This interpretation of Fig. 2 is not much altered by varying the parameters relevant with 
f. First cross sections of spallation processes are inversely proportional to mf, and then the 
timescale of each process linearly increases as mf increases. Thus, even when mf is larger than 
mf = 350GeV by up to a factor of ten, the region of 6 Li over-production scarcely changes. Next 
we point out that our result depend only mildly on the left-right mixing of the stau. Indeed, 



cross section of the 4 He spallation is proportional to sin 2 6 T . Its order of magnitude will not 
change as long as the right-handed component is significant. 

4. Summary 

Long-lived charged massive particles provides some exotic nuclear reactions in the big bang 
nucleosynthesis. So it is important for understanding the property of long-lived charged massive 
particles to understand what type of exotic nuclear induce over-production (-destruction). Newly 
included in the present work is the spallation of the 4 He in the stau- 4 He bound state given in 
Eq.(l). This process is only present in the model which predicts the long-lived charged particles 
due to the phase space suppression with the weakly interacting daughter particle. 

We calculated the rate of the helium-4 spallation processes analytically, and compared it with 
that of catalyzed fusion. We found that the spallation of 4 He nuclei dominate over the catalyzed 
fusion as long as the phase space of the spallation processes are open and hence the property of 
long lived stau is constrained from avoiding the overproduction of a deuteron and/or a triton. 
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